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Superfluidity in atomic Fermi gases with population imbalance has recently become an exciting research 
focus. There is considerable disagreement in the literature about the appropriate stability conditions for states 
in the phase diagram throughout the BCS to Bose-Einstein condensation (BEC) crossover. Here we discuss 
these stability conditions for homogeneous polarized superfluid phases, and compare with recent alternative 
proposals. The requirement of a positive second order partial derivative of the thermodynamic potential with 
respect to the fermionic excitation gap A (at fixed chemical potentials) is demonstrated to be equivalent to the 
positive definiteness of the particle number susceptibility matrix. In addition, we show the positivity of the 
effective pair mass constitutes another nontrivial stability condition. These conditions determine the stability of 
the system towards phase separation of one form or another. We also study systematically the effects of finite 
temperature and the related pseudogap on the phase diagrams defined by our stability conditions. 
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I. INTRODUCTION 

Superfluidity in atomic Fermigases has become an impor- 
tant research arena QJ 00015100 Hi for both the condensed 
matter and atomic, molecular and optical physics communi- 
ties. In these ultracold gases, via a Feshbach resonance, one 
can tune the pairing interaction strength continuously from 
very weak in the BCS limit to very strong in the Bose-Einstein 
condensation (BEC) limit. Adding greatly to the excitement 
has been a recent emphasis on experimental studies of super- 
fluidity in the presence of a population imbalance 
between the two fermion spin components. This system has 
important consequences for other subfields of physics includ- 
ing nuclear physics and high density QCD ILL?! fl3l. 11411 . From 
a theoretical standpoint this problem is particularly rich and 
at the same time complex lITafTrllr^lmlT^IMEnElEl . 
Simple mean-field calculations 1 24] show that, unlike in the 
equal spin mixture case, a homogeneous superfluid state is not 
always a stable ground state at zero temperature (T). There 
are a number of different ground states to consider. More- 
over, there has been considerable controversy in the literature 
l!2lll24ll25Ll26il about the precise nature of the stability condi- 
tions, associated with one mean-field ground state or another. 

It, thus, becomes particularly important to characterize 
and study the stability conditions associated with polarized 
superfluids. Moreover, since there is a natural extension 
HlHlilllHHii of the standard ground states to finite tem- 
peratures T, in this paper, we derive these stability conditions 
for a homogeneous interacting Fermi gas superfluid with pop- 
ulation imbalance throughout the entire BCS-BEC crossover, 
and at arbitrary T. We will focus our attention on a super- 
fluid composed of a condensate of pairs which has a zero total 
center-of-mass momentum. Therefore, we will not discuss the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state BIB , which 
allows condensation of pairs at finite momenta. It appears that 
the FFLO state with pairing at only one value of momentum 
qo is stable only in a very limited phase space 13211 . Multiple 
plane wave FFLO states are much more complicated and are 
currently under investigation. 



We set up the central issues of this paper by summarizing 
our previously obtained 1 24] results on the zero temperature 
phase diagram in the p-l/kpa plane plotted in Fig.[T] Here 
p = 5N/N is the polarization, kp is the Fermi wavevectorof a 
noninteracting Fermi gas of the same number density without 
a population imbalance, a is the two-body s-wave scattering 
length, N = + N L and 5N = 7V T - 7V X > are the total 
number and number difference, respectively. In the BCS limit, 
our result for the boundary separating the stable polarized nor- 
mal Fermi gas and the unstable Sarma phase agrees with that 
in Ref. 1 25], but differs quantitatively from that in Ref. 1 26] . 
On the BEC side, our result agrees with that of Ref. 12611 but 
not that of Ref. 1 25 ] . Since there is so much controversy in 
the literature, this paper, then, addresses a particularly timely 
issue. 

The differences between the various theoretical proposals 
for the phase diagram stem from different conclusions con- 
cerning the stability of the various phases. There seems to be 
general agreement I2lll24l.l26ll about the two generic stability 
requirements which were first articulated in Ref. 1 25] where 
it was argued that the number susceptibility matrix must have 
positive eigenvalues, and that the superfluid density must also 
be positive. We will show in detail here that in contrast to 
Ref. 1 26], that the positivity of eigenvalues of the number sus- 
ceptibility matrix is equivalent to the positivity of the 2nd or- 
der partial derivative J^J, where Q is the thermodynamic po- 
tential, and A the fermionic excitation gap. The same obser- 
vation was very recently made in Ref. |21]. In comparison 
to these two equivalent conditions, we find that the positivity 
requirement on the superfluid density is much less st ring ent. 
This differs from Ref. [21] but agrees with Refs. 1251 12611 . In 
addition, we find that the requirement that the pair mass be 
positive constitutes another nontrivial stability condition, and 
may be more stringent than the positivity of at low tem- 
perature and low imbalance in the BCS regime. 

Our theoretical formalism is described in Ref. 1 24], so we 
will not repeat the details here. The system is composed of 
two spin components of number N-\ and JVj, and of chemical 
potential /if and /ij, respectively. They interact via the attrac- 
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FIG. 1: Zero temperature phase diagram of a homogeneous Fermi 
gas as a function of pairing interaction characterized by l/fe^a and 
the polarization p. On the left side of the (black) T C MF = curve, 
the system is a stable polarized Fermi gas (labeled by "FG"). On the 
right side of the (green) d 2 Q/dA 2 — curve, the system is a stable 
polarized superfluid at low T, labeled by "SF". At T — 0, the entire 
shaded region between these two curves is unstable against phase 
separation, and a stable polarized superfluid can exist only in the 
BEC regime. The (red) n s = curve and the (blue) T c = = 1 /M* 
curves are completely within the unstable regime. Here kp is the 
Fermi wavevector of a noninteracting two component Fermi gas of 
the same density n without a population imbalance. All energies are 
measured in units of Ef = k%/2m. Throughout this paper, we use 
fco = 80k F . 



tive interaction E/k.k' = U <^k<^k with U < 0. Here we use a 
Gaussian cutoff tp^ = er k l 2k ° with fco sufficiently large (as 
appropriate for a short range potential). The cutoff t^k can be 
used to approximate a potential, which includes but also gen- 
eralizes the well studied contact potential. Here fco is given by 
the inverse range of interaction, and we take fco = 80fcF in our 
numeric calculations throughout this paper. This interaction 
is related to 1/kpa via m/Aira = 1/U + ^ k <^ k /2ek, where 
£k = k 2 /2m. The full fermion Green's function is given by 
G] t i(K) = u k / (iuj n ±h— i?k) +«k/ ( iuj n ±h + E}k), where 

= v^k + A2 ' & = £k - H, M = (MT + Ml)/ 2 ' £_J= 
(MT - W)/ 2 and U L v l = C 1 ± Ck/-Bk)/2. As in Ref. [24], 
we set the volume V" = 1, H = ks = 1, and if = (iw n , k), 
Q = (iQ m ,q), J]k = TY.„ V,,, etc, where w n (O m ) is 



odd (even) Matsubara frequencies 1 33 ] . Our Green's function 
is equivalent to presuming the BCS form for the self-energy, 
E a (K) = -A 2 G°(-K)tpl, where a =t,| and a = -a. 
Here the bare Green's function Gq\{K) — iui n — ^k,o-- With 
these definitions we introduce the important pair susceptibility 
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At finite T, we have four equations in the presence of pop- 
ulation imbalance. They are the gap (A) equation, two num- 
ber equations, and the pseudogap (A p9 ) equation. At zero 
T (where A pg = 0), our theory yields three equations and 
reduces to the standard one-channel results in the literature. 
More generally A 2 (T) = A 2 SC (T) + A^(T), where A sc is 
the superfluid order parameter. Except for the number differ- 
ence equation, all equations can be written in the same form as 
in the case with equal spin population [28], provided one re- 
places the Fermi distribution function f(x) and its derivative 
f'(x) with the averages f(x) = [f(x + h) + f(x — h)}/2 and 
f'{x). At T ^ 0, we must include the effects of noncondensed 
pairs via the pseudogap term A pg . These noncondensed pairs 
have been ignored in previous work ll2lll34ll35ll . and are found 
to be extremely important here. They have an effective pair 
mass M* which must necessarily be positive, thereby, adding 
another condition to the requirements for phase stability. 

At general T, our self consistent equations are 

= l + U X (0) = l + Uj2 1 ~* f E {Ek) vl, (2) 
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pn = J2[f(E k -h)-f(E k + h)] 



Q#0 q 
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which enters throughout the present paper. We thus have 

J2 K Z a (K)G a (K) = -A 2 X (0), 



Here x(0) is the pair susceptibility at Q = 0, and f2 q = 
q 2 /2M* is the pair dispersion. From Eq. (I3bt . we have h > A 
at T = for any p > 0. Solving Eqs. Q-Q with A = at 
T = 0, we can obtain the boundary between the polarized 
normal Fermi gas phase and the unstable phase in Fig.Q 



II. THERMODYNAMIC POTENTIAL fi 

We first discuss the gap [Eq. @] and number equations 
[Eqs. yJ], which govern the fermionic degrees of freedom, 
and only later address pseudogap effects which appear at fi- 
nite T through Eq. @. The gap and number equations can be 
obtained from the fermionic part of the thermodynamic poten- 
tial fi(A, /j,, h) via partial derivatives with respect to A, \x, and 
h, respectively. It should be noted that only the gap equation 
corresponds to a vanishing first order derivative. 

Instead of writing down the thermodynamic potential from 
the known quasiparticle energy spectra 1 16], we calculate it 
via the energy and entropy. We have 

K = E — /if iVj - n i N i = E - fiN - hSN, (5a) 
fi = E - TS - fiN - hSN = K - TS , (5b) 
F = E -TS = Q + fiN + hSN , (5c) 
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Additional second order partial derivatives of interest are 

d 2 n 



A, fj. 



The energy for each spin is given by 
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G(r(K) . 
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Using the gap equation (|2}, one can show that 

E = J2E« =E e k[G T (/O + G i (K)]-A 2 x (0) 
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K = E - ^ + 2^k/(i?k)] -hSN- — . (9) = _^ = 2 £ jf (£ k) 



The quasiparticle excitation energies are now given by 
±i?k + c/i, where cr = ±1 for spin down and up, respectively. 
The entropy contains contributions from both spin species. 



S = - E Et/(-^ + ^)ln/(-i?k + ^) 

CT=±1 k 

+ f {E k + ah) In f(E k + ah)} . 
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III. STABILITY CONDITIONS 
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Finally, we obtain 
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(ii) 



where we have used the relationship dE^/dA = A(p k /i? k . It 
is evident that when the gap equation is satisfied, dVl /dA — 
0. Alternatively, one easily finds dF/dA — at the same 
time. It is interesting to note that in Eq. ( Bel the term /j,N 
cancels that from E, whereas hSN cancels that in TS. It is 
straightforward to verify Eqs. (|6} are equivalent to the usual 
form of the number and number difference equations d3al and 
d3bl . Using Eqs. ©, 

When the gap equation is satisfied, we readily obtain 
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A. Stability of the polarized normal phase on the BCS side 

The condition dil j dA = admits a trivial solution A = 0. 
Referring now to Fig. [2 we see that on the left side of the 
boundary T~* = between the A = and the paired phase, 
the quantity %(0) is a function of \i and h only. Here we have 
U~ x + x(0) < 0, independent of 1/kpa. It is easy to show 
that the expression inside the square brackets of the second 
line in Eq. (I12a> vanishes when _E k < h. Therefore, when 
A = but 1 + Ux(0) ^ 0, we always have 



d 2 il 
dA 2 



= -2 



1 

U 



x(o) 



> 0. 



(13) 



We will argue later that the positivity of is associated with 
stability against phase separation. From this and the above 
inequality, we can conclude that the homogeneous polarized 
Fermi gas phase in Fig.^is stable with respect to phase sepa- 
ration. 



d 2 F 

i. Second order total derivative — > - Stability of 
pairing versus normal Fermi gas for fixed N and SN 



Next we investigate the total derivative 



d 2 F 



dA 



2 ' 



for a fixed 



particle number system. This is equivalent to confining the 
equations in the hyperplane defined by the two number equa- 
tions in the parameter space spanned by A,fi and h. The 
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constraints imposed by the number equations Q exclude the 
possibility of phase separation or any spatial variation of the 
order parameter. We will show that d 2 i^/dA 2 > provides 
a stability condition for superfluidity vis a vis a polarized nor- 
mal Fermi gas phase. This is equivalent to the statement that 
the free energy reaches a minimum when a nontrivial solution 
(A ^ 0) of our equation set is found. 
From the number equations, we have 



dN _ dN dN dN dh 
dA _ dA + dA + ~dhdA ' 
dSN dSN dSN d^ dSN dh 



so that 
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FIG. 2: Behavior of d 2 F/dA 2 at T = as a function of polariza- 
tion p for various interaction strength 1/kpa, from BCS to BEC, as 
labeled in the figure. From top to bottom, 1/kpa decreases monoton- 
ically. d 2 _F/dA 2 remains positive on the right side of the T^ IF — 
boundary in Fig.Q and vanishes on this boundary. 



Using the chain rule, 
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Number susceptibility and stability of homogeneous 
polarized superfluid against phase separation 



finally, we obtain 
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where we have used Eq. ( 1171 to simplify this expression. 

We plot the total derivative d 2 -F/dA 2 at T = as a func- 
tion of polarization in Fig. [2 This figure demonstrates that for 
all regions on the right side of the T C MF = line in Fig. Q 
(where the solution of the gap equation with fixed particle 
number exists), d 2 F/dA 2 is always positive. It vanishes pre- 
cisely at this boundary, that is, on the line T C MF = 0, which 
separates the paired phase and the polarized normal Fermi gas 
phase. For definiteness, at unitarity, d 2 F/dA 2 vanishes at 
p = 0.72. Above this polarization the system is stable as a 
normal, unpaired Fermi gas. Extrapolation of d 2 F/dA 2 to 
the normal Fermi gas side of this line would lead to a negative 
value, indicating that the paired state is not stable there. 

In summary, when a solution of the self-consistent equa- 
tions (|2ji-@ can be found, a paired phase is expected to be 
more stable than the unpaired Fermi gas phase, both at zero 
and finite T. This result agrees with what we obtained in the 
previous subsection llll Al but addresses this issue from a dif- 
ferent perspective. Here we consider the stability of the paired 
phase whereas in Section ITlI Al we addressed that of the un- 
paired Fermi gas. We can turn this around to conclude that the 
positivity of 'd 2 F ' j 'dA 2 does not provide an extra constraint on 
the stability of the polarized superfluid phase. 



In this section, we consider the number susceptibility with 
respect to variation of the chemical potentials in equilibrium. 
It is expected on physical grounds 1 25] that both eigenvalues 
of the matrix dN a /d/j, a ' must be positive in order for the sys- 
tem to be stable, against phase separation. It should be obvi- 
ous that when A is held fixed, the eigenvalues of the matrix 
dNfj/d/jLa' are always positive. However, when A is treated 
as an implicit function of // and h as defined by the gap equa- 
tion (0, the eigenvalues may change sign. 

For simplicity in notation, we use partial derivative d/dx 
when A, n and h are all treated as independent variables, as 
in the previous section. When A is regarded as a function of 
/i and h, we define 
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dA d 
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(x = n,h). 
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With the relationships N = N±, SN = — N±, and 
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Here a = a - x = m\ -i 



is a unitary transformation 
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matrix. Therefore, the eigenvalues of the matrix 
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are also required to be positive. 

From the gap equation, we readily obtain 
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It should be noted that this expression is valid only when 
the gap equation is satisfied, i.e., on the right side of the 
T C MF = boundary in Fig. □ On the left side, we have the 
inequality of Eq. (1131 . 



Finally, we have 



DN 


d 2 n 
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Therefore, the eigenvalues are given by 

Tr(M 2x2 ; 



A+ = 



± v/Tr(M 2x2 ) 2 



4det(M 2x2 ) 



(25) 



where Tr(M 2x2 ) 



DJV 



DSN 
Dh 



and det(M 2 



x2) 



DN DSN 
D/i Dh 



are t ' le trace an d determinant, respectively. 



Since appears in the denominator of the expressions in 

Eqs. J24L each element of M 2x2 will change sign when 
approaches zero and changes sign. Under this circumstance, 
[det(M 2x2 )] _1 approaches zero and changes sign. Our nu- 
merics shows that only one of the two eigenvalues in Eq. (125 \ 
changes sign. This roughly corresponds to ^jf. Therefore, 
the stability condition that the eigenvalues of the number sus- 
ceptibility matrix be positive is equivalent to 



d 2 n 

SA 2 



> , when 1 + Ux(0) = . 



(26) 



This condition has been argued in the literature to corre- 
spond to the condition for phase separation. Indeed, this con- 
dition says when fi and h are held fixed, a stable homoge- 
neous superfluid solution should always minimize the ther- 
modynamic potential. Otherwise, the system tends to phase 
separate into a region with smaller A and another region with 
larger A, both of which give lower SI. In such a case, the num- 
ber density in each region is not fixed. Such a phase separation 
is generally believed to be realized by a one-component Fermi 



gas physically adjacent to an unpolarized superfluid. We must 
have A > ft, in the unpolarized superfluid region. If we take 
into account the possibility of condensation at qo ^ 0, this 
form of microscopic phase separation may not be the only way 
to address the instability. One may have a more homogeneous 
state as well such as found in FFLO-like phases. 

The behavior of <9 2 S1/<9A 2 , the inverse determinant 
1/ det(M 2x2 ) = 1/(A+A_) and the eigenvalues of the num- 
ber susceptibility matrix M 2x2 at T = are all plotted in 
Fig. [3] as a function of polarization p for various interaction 
strengths 1/kpa. From Eqs. J24> . it is evident that all ele- 
ments of the matrix M 2x2 change sign where <9 2 S1/<9A 2 = 
0. Nevertheless, we notice that one of the two eigenvalues, 
EVi = A+, is always positive. This roughly corresponds 
to DN/D/j,. However, the second eigenvalue, EV2 = A_, 
does change sign for 1/kpa — 1, 1.5, and 2, exactly where 
d 2 fl/dA 2 = 0. Moreover, the fact that d 2 fl/dA 2 appears in 
the denominator in Eqs. (1241 is manifested by the jump in X± 
upon the sign change. 

If one extrapolates d 2 il/dA 2 in the upper left panel, for 
1/kpa < 0.3, one would see that this partial derivative be- 
comes positive above the T^ IF = line in Fig.^ This recon- 
firms that the polarized normal Fermi gas phase in Fig. ^ is 
stable. One can see thus that the T^' IF = provides another 
d 2 fl/dA 2 = (at A — > 0) line, but this time on the BCS side 
of the resonance. Note the contrasting behavior of d 2 £l/dA 2 
and d 2 F/dA 2 across the T c Mi? = line. 

At this stage it is useful to compare with the literature. It 
is argued in Ref. 1 26] that the positivity of eigenvalues of the 
matrix Af 2x2 is a necessary but not sufficient condition for 
stability. Moreover, in this previous work, the positivity of 
eigenvalues of the number susceptibility matrix was found to 
be different from the condition contained in Eq. i26\ . By con- 
trast, here we find that these two conditions are indeed equiv- 
alent. The same observation was made in Ref. 12111 . although 
only the matrix element DSN/Dh was addressed. 

In other recent work a numerical optimization pro- 

cedure was invoked in which the entire landscape of il as a 
function of A for fixed fi and h was used to find the most sta- 
ble phase. It was argued that d 2 il/dA 2 = often gave multi- 
ple solutions, sometimes corresponding to a local maximum. 
In such cases, the A = solution was sometimes claimed to 
be more stable, since it minimizes globally. 

These results should be contrasted with the present calcu- 
lations. Except for the trivial A = solution, we find there 
is at most one solution to our set of equations at either zero or 
finite T, We argue that the optimization method is in principle 
valid, but it requires that the system be in chemical equilib- 
rium with an infinitely large particle reservoir so that \x and 
h are unchanged before and after a possible phase separation. 
Because the atomic Fermi gases constitute a finite system, p, 
and h will be different before and after phase separation; they, 
thus, no longer satisfy the same particle number constraint. In 
this way we conclude that one cannot use a single SI (A) curve 
to find the solution for the stable phase for a given polarization 
at fixed particle number. 
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FIG. 3: Behavior of d 2 tt/d A 2 , the inverse determinant 1/ det(M2 X 2) and the eigenvalues of the number susceptibility matrix M'2 X 2 at T = 
as a function of polarization p for various interaction strength l/fefci, from unitarity to BEC, as labeled in the figure. The first eigenvalue, 
EVi = A+, remains positive, roughly corresponding to DN/Dfi. The sign change of d 2 Q./dA 2 and the second eigenvalue, EV2 = A_, and 
thus 1/ det(M2x2) occurs simultaneously, for 1/kFO, = 1, 1.5 and 2 on the BEC side of the Feshbach resonance, along the d 2 fl/d A 2 = 
phase boundary in Fig.Q 



D. Superfluid density 

For a superfluid state to be stable, another obvious condi- 
tion ll25ll is that the superfluid density must be positive. Us- 
ing linear response theory, one can calculate the superfluid 
density via the response of the system to an external (ficti- 
tious) vector potential, as if it were a charged superconductor. 
In this way one computes the the current-current correlation 
functions compatible with our T-matrix approximation 12811 . 
as shown in Ref. (36ll . Here, however, we need to keep track 
of the population imbalance. Without showing the details, our 
superfluid density is given by 



than that given by Eq. i26\ . corresponding to the positive def- 
initeness of the number susceptibility matrix. This is different 
from Ref. 

It should be noted that at finite T > T c , n s vanishes identi- 
cally. However. n^ IF may remain finite as long as the mean- 
field equations and Q are satisfied. In this case, nf tF may 
change sign. 

Our numerical result shows that the n s = line shifts to the 
upper left in the p—l/kpa phase diagram as T increases from 
zero, and disappears completely slightly below T jTp < 0.2. 
Based on Fig. we conclude that n s remains non-negative 
for l/fc^a > 0.5. 
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^SC 

A 2 



> 0. 



(27) 



where n^ IF is an artificial construct corresponding to the su- 
perfluid density in a BCS-like strictly mean-field treatment. 
The boundary n s = at T = is given by the red line 
in Fig. ^ This line is completely in the fermionic regime 
(p > 0) and within the phase boundary set by T^ IF = 
and d 2 fl/dA 2 — 0. The superfluid density is positive on 
the right side of this line, and negative otherwise. We may 
summarize. The positivity of n s is a much weaker constraint 



E. Positivity of the effective pair mass 

The effective pair mass can be obtained by Taylor expand- 
ing the important pair susceptibility x(Q) m me form 



x(Q) - x(o) = z m m - 



2AP 



(28) 



Most theoretical work in the literature on polarized Fermi 
gases is based on the ground state El, or at most on a BCS- 
like mean-field treatment at finite T I2lll35ll . which does not 
include noncondensed pairs. Here, by contrast, we include 
these noncondensed pairs (A 2 g ^ 0) or pseudogap effects 



II27II which are reflected in the pair propagator (or T-matrix) 
and its associated dispersion. As a result, we require that the 
pair mass be positive. This condition appears as another line 
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FIG. 4: Phase diagram in the T — 1/kpa plane for representative values of p = 0.1, 0.15, 0.2, and 0.5 (as labeled). Shown are T c (blue line), 
T C MF (black), and the instability boundaries defined by d 2 Q./dA 2 — (green), and l/M* = (red line), respectively. The yellow shaded 
area represents stable superfluid (labeled by "SF"), the dotted region is unstable. The white open space on the right of the T C MF line represents 
a stable pseudogap phase ("PG"), whereas on its left lies the stable unpaired Fermi gas state ("FG"). At p m 0.14, the superfluid region splits 
into two, and the small closed region shrinks with increasing p, and disappears at p ~ 0.18. Forp = 0.1 and 0.15, the two instability lines 
intersect with each other, whereas for p — 0.2 and 0.5, l/M* = is completely inside the unstable region defined by d 2 fl/dA 2 < 0. The 
(segment of each) stability line appear as a dashed line when it is inside an unstable region defined by the other stability condition. We also 
show the line on which n s = (cyan dashed curve) but for p = 0.5 only since it is always inside the unstable region. 



in the phase diagram of Fig. G] At T = 0, this would corre- 
spond to the blue T c = line, which, at this temperature, is 
completely inside the unstable regime and does not provide an 
extra constraint. However, at finite T, this line moves consid- 
erably to the BCS side, especially at low p. As will be shown 
below, at low T and low p, this imposes a more stringent con- 
straint than Eq. (I26> in the BCS regime. 



IV. TEMPERATURE DEPENDENCE OF THE PHASE 
DIAGRAM 

Since superfluidity is essentially a finite temperature phe- 
nomenon, it is important to assess the stability of the various 
superfluid phases at T ^ 0. This is where experiments are 
performed. Our theory allows us to calculate the superfluid 
transition temperature T c in the presence of the pseudogap ef- 
fects associated with incoherent finite center-of-mass momen- 
tum pair excitations. We will see that quite systematically, su- 
perfluid phases which are unstable at zero temperature, may 
become stable at finite T. As discussed in Ref. 1 24] these in- 



termediate temperature superfluids will behave somewhat dif- 
ferently depending on whether they occur to the right or the 
left of the T c = quantum critical point line in Fig. [2 

In Fig.0]we present a slice of the finite temperature phase 
diagram, plotted as characteristic temperature versus 1/kpa, 
at four representative polarizations, p = 0.1, 0.15, 0.2 and 0.5. 
The lines in the figures correspond to T c and the phase bound- 
aries defined by the instability conditions, including where 
pairs are no longer found to be stable (via l/M* < 0). The 
shaded regions indicate where stable superfluidity is found. 
The dotted regions indicate where either the phase separation 
stability criterion (given by > 0) is violated, or the effec- 
tive pair mass becomes negative. It is important to note that 
on the BCS side of resonance (depending on the polarization) 
there may be two values of T c for each value of 1/kpa. The 
meaning of these two T c 's was discussed in earlier work 1 24]. 
This case applies to the unitary regime, as should be seen in 
the figure. Below the lower T c we find no solution to our set of 
four equations (|2j-©. In this way the system cannot support 
superfluidity. More concretely, the effective pair chemical po- 
tential become negative below the lower T c . (At even lower 



8 



0.8 



0.6 



0.4 



0.2 



d(fcpa)/d7" c = 
3 2 0_/3A 2 = at T c 

<y'(.i r).\* o /:::•: 

tic / 1 ° ' ." 

r ^ =0 / : 11 

/. o /:•!►?:•: : : :f v 

. j — ► 
■IS 

'da/dT^O ^ 



0.3 



1/^ = 0.2 



-3-2-10123 
Mk F a 

FIG. 5: Phase diagram in the p — 1/kpa plane, showing where 
intermediate temperature superfluidity (shaded region) exists. The 
T C MF = 0, T c = 0, and zero temperature d 2 Q/dA 2 = lines are 
the same as in Fig.Q The red curves represent where d 2 Q/dA 2 — 
at corresponding T c . The line defined by A(l/kFa)/AT c = is 
given by the turning points (p, 1/kpa) where 1/kpa reaches a local 
extremum as a function of T c . At T — 0, the entire region between 
the T C MF = and the d 2 Q,/dA 2 = lines is unstable against phase 
separation. However, a stable polarized, intermediate temperature 
superfluid phase exists for the yellow shaded region enclosed by the 
(orange) A(kpa)/AT C = curve and the p = axis, or by the (red) 
d 2 Q(T = T c )/dA 2 = and (green) d 2 Q./dA 2 = curves. Super- 
fluidity in the dotted region is unstable at any temperature. 




T, the effective pair mass becomes negative). Above the up- 
per T c , superfluidity is destroyed in the usual way by finite 
temperature effects and the system becomes a normal Fermi 
gas, which is far from ideal. Here there are still strong pair- 
ing correlations giving rise to a pseudogap in the fermionic 
excitation spectrum. 

By contrast as we move towards the BEC regime we find 
only one T c , but the stability of the superfluid at low T is 
cut off because of the negativity of Finally, sufficiently 
deep into the BEC regime, a stable superfluid persists for all 
temperatures (below T c ) including T = 0. 

These results are summarized in the form of a finite tem- 
perature phase diagram in Fig. [5] The shaded region indicates 
where intermediate temperature superfluidity occurs, while 
the dotted region corresponds to an unstable superfluid, pre- 
sumably a LOFF-like or heterogeneous phase. Because there 
is a regime where two T' c s are evident, as seen in Fig.0 we 
indicate the upper boundary of the shaded region by a "turning 
point condition": d(l/kpa)/<lT c = 0. 

The condition given by Eq. (I26> at T c is represented by the 
red curve in Fig. [5] Note that T c self-consistently along this 
line. 

The next two figures correspond to regions on either side of 
the quantum critical point T c = in Figs.ffland [5] The region 
to the right of this line ("Regions IIB and IID" in Ref. 12411 ') 
is further discussed in Fig.|6] It corresponds to the BEC side 
of resonance above 1/fepct = 0.2. On the left of this line 



FIG. 6: Behavior of T c and the instability phase boundary on the 
BEC side. T uns t a Ue is defined by d 2 fl/dA 2 = in the T-p plane, 
for 1/&FCE = 0.2, 0.3, and 1, from top to bottom. For each value 
of 1/kFa, the phase diagram is composed of four different phases, 
separated by the solid lines: stable superfluid (SF), unstable super- 
fluid (uSF), stable pseudogap (PG), and unstable pseudogap (uPG) 
phases, as labeled in the figure. The yellow shaded region, on the 
left of the (black) T c and (red) T unsta ti e curves, e.g., the shaded area 
for 1/kFa = 0.2, is a stable polarized superfluid. The superfluid so- 
lution below these two curves is unstable. Above these two curves, 
there exists a stable pseudogap phase. The unstable phases will dis- 
appear when 1/kpa > 2.3. The pair mass is always positive for 
ljkpa > 0.2. We also show for 1/kpa = 0.2 the n 3 — line, 
which is completely within the unstable region defined by Eq. 1261 . 
n s does not change sign at any T and p for 1/kpa > 0.5. 

("Regions IIA and IIC" in Ref. Q) which includes the BCS 
side of resonance as well as the unitary phase we summarize 
our findings in Fig.Qbelow. 

In Fig. [6] we show how T c and T unsta bie evolve with p and 
1/kpa on the BEC side, where T unsta Ue is the temperature 
where ^| = 0. We do not show the deep BEC case where 
the superfluid phase is always stable. Here we plot these two 
families of curves for 1/kpa = 0.2, 0.3, and 1.0, from top to 
bottom. The upper and lower curves in each panel are for T c 
and for T unsta bie, respectively. Each phase diagram is com- 
posed of four different phases. Inside the shaded area, above 
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FIG. 7: Phase diagram on the T-p plane on the BCS side of the 
T c = line in Fig.Q for l/k F a = 0.1, 0, and -0.5, from top to 
bottom. In each panel, the superfluid region is shaded in yellow, on 
the left side of the T c curve. Except for 1/kpa = 0.1, the (red) 
T uns tabie line does not intersect with the T c curve. There is a small 
unstable superfluid phase (uSF) for l/kwa = 0.1, as labeled in the 
top panel. For X/kpa = —0.5, thepseudogap (PG) is present only in 
a narrow temperature range since here the pair formation temperature 
T* (approximately given by T C MF ) is not much higher than the upper 
T c . At higher T, the phase is a polarized normal Fermi gas. Slightly 
below the lower T c , the pair mass becomes negative (below the green 
1/M* = line), implying that the paired phase in the mean-field 
treatment is unstable in this region. The AT > requirement is a 
tighter constraint than Eq. 1261 in the BCS regime at low T and low 
P- 



T U nstabie but below T c , there is a stable intermediate tem- 
perature superfluid. Below these two curves, the superfluid 
phase is unstable and this is denoted by "uSF". To the right 
of the T c curve we have a stable pseudogap phase above the 
Tunstabie line, labeled by "PG". Below the T unsta Ue line we 
have an unstable pseudogap phase labeled "uPG". When there 
is a PG phase (stable or unstable) the pairing onset tempera- 
ture T* (approximately given by T C MF , not shown) is higher 
than T c . As 1/kpa increases, the intersection point between 
Tunstabie and the T = axis moves to the right, and the 
area of the stable superfluid phase grows, until 1/kpa « 2.3, 



where the area of the unstable phase shrinks to zero. The 
case 1/kpa = 0.2 is unusual; the T c curve bends to the 
left at low temperature. We also show the n s = line for 
1/kpa = 0.2. It is completely within the unstable regime, 
and it disappears for 1/kpa > 0.5. The pair mass is always 
positive for 1/kpa > 0.2. 

It should be noted that at a given T below the point where 
T c = Tunstabie, there is a critical polarization p c , above which 
the superfluid becomes unstable against phase separation (or 
an FFLO-like pairing state). From Fig. [6] it is evident that p c 
increases with T. This trend is consistent with recent experi- 
mental results reported by the Rice group 13711 - which report 
that the critical polarization for phase separation decreases 
with decreasing T. In the present work, of course, we have 
not incorporated the effects of the trap potential. These were 
discussed elsewhere 1 30 ] . 

We now address Fig. which presents similar plots in the 
T-p plane for fixed values of 1/kpa = 0.1, 0, and -0.5. This 
corresponds to the left side of the T c = curve in Fig. ^ 
Generally, the T uns tabie curve lies below the corresponding 
lower T c curve. The shaded region on the left of the T c curve 
in each of the three panels is a stable superfluid phase. For 
1/kpa = 0.1, however, the superfluid phase has a tiny region 
of instability enclosed by the T c and T un stabie curves. As in 
Fig. |6j the normal region is split into a stable and unstable 
phase, where pseudogap effects persist except at high T. In 
all cases, the paired phase below T un stabie is unstable. The 
effective pair mass becomes negative below the line defined 
by 1/M* = 0, which intersects the T un stabie line. Here the 
positivity of the pair mass provides a more stringent constraint 
on stability than Eq. J26I . Finally, just as in Fig.|S] one can 
find a critical polarization p c for these values of 1/kpa as 
well, and p c increases with T. 

It should be noted that Figs. [5] and evolve smoothly into 
each other. As 1/kpa changes continuously from 0.2 to 0.1, 
the T c curve at low T in the top panel in Figs. [6] will bend 
further to the left, and eventually touch the p = axis and 
change into the top panel of Fig.Q 

We end with Fig. [8] This shows how the various phase 
boundaries first presented in Fig. Revolve with temperature. 
The three panels, from top to bottom, correspond to three 
different temperatures T/Tp = 0.05, 0.1, and 0.2, respec- 
tively. The solid lines split the phase diagram into a few dis- 
tinct phases. The points at which T" c MF is equal to the given 
temperature of each figure lead to a line which separates the 
unpaired Fermi gas (FG) phase from the paired phases. Sim- 
ilarly, the T c — T line separates the superfluid phase (be- 
low the line) from the pseudogap or normal phases (above 
the line). Finally, the lines associated with <9 2 f2/<9A 2 = 
and 1/M* = separate the stable superfluid and pseudogap 
phases (below this line) from unstable phases. These two lines 
intersect with each other, and when this happens, one will ap- 
pear partly inside the unstable region defined by the other, as 
indicated by the red and green dashed lines. Thus the sta- 
ble superfluid phase lies in the yellow shaded region. The 
n^ IF = line (not shown) always appears completely within 
the unstable phases, and therefore, does not provide a separate 
phase boundary, as in the T — case. (Note here we use n^ IF 
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FIG. 8: Evolution of the phase diagram in the p — 1/kpa plane with 
temperature. From top to bottom, T/Tf = 0.05, 0.1, and 0.2, re- 
spectively. The phase diagrams are split into different phases by the 
solid curves, as labeled in the figure. Except the n = line (shown 
for T/Tf = 0.2 only) stay relatively unchanged, all other phase 
boundaries moves with T. The mean-field "phase transition" line, 
defined by T C MF = T, moves to the right. The (blue) T c = T line 
separates the superfluid and the pseudogap (or normal) phase, and the 
(green) d 2 fi/ dA 2 — line, in conjunction with the (red) 1 /M* = 
line, separates the stable and unstable phases. Stable superfluid exists 
in the yellow shaded region, the dotted region is either an unstable 
superfluid (labeled by "uSF") or unstable pseudogap ("uPG") phase. 
The white space on the right of T C MF = T line is a stable pseudo- 
gap ("PG") phase. The 1/M* > represents a stronger constraint 
than Eq. 1261 only at low T and low p (e.g., in the top panel). The 
1/M* = line moves to the upper right and disappears at slightly 
below T — 0.2. The sign change of rig (not shown) always occurs 
within the unstable regions. We also show /i = which separates 
fermionic from bosonic regimes. 



instead of n s because n s — outside the superfluid regions.) 

As temperature increases, the line associated with F = 
T moves to the right. At the same time both the 1 /M* = 
curve and the n^ 1 = curve shrink and move to the upper 
left, ultimately disappearing slightly below T/Tp = 0.2. The 
fi = boundary between fermionic and bosonic regimes is 
rather insensitive to T, and so is indicated just for the case of 
the highest temperature. 

The evolution of the T c = T and d 2 Q/dA 2 — curves is 
rather interesting. For low polarization p, both curves move 
to the left into the BCS regime. This reflects the fact that 
temperature stabilizes the paired and superfluid phase. This 
is consistent with our observation of intermediate temperature 
superfluidity in the presence of population imbalance. The 
T c = T curve on the BEC side is continuously suppressed 
by raising T. This is simply due to fact that T c decreases 
with p for given 1/kpa in this regime. In contrast with the 
T = case, we have a stable pseudogap phase appearing at 
finite T. This corresponds to the white region on the right of 
the T™ F = T line, below the d 2 n/dA 2 = and 1/M* = 
lines, and above the T c = T line. At T/T F = 0.2, the 
d 2 fl/dA 2 — and the T c = T lines do not intersect each 
other, so that the entire superfluid phase is stable. 



V. CONCLUSION 

In summary, in this paper we have studied the stability con- 
ditions for polarized fermionic superfluids in considerable de- 
tail. We have restricted our attention to zero momentum con- 
densate pairs, thereby, excluding FFLO-like phases. We find 
that the positive definiteness of the number susceptibility ma- 
trix is equivalent to the positivity of (d 2 il/dA 2 ) fl provided 
the gap equation (|2ji is satisfied. In addition, we find that the 
positivity requirement of the effective pair mass constitutes 
another nontrivial stability condition. At low temperature and 
low population imbalance in the fermionic regime, this latter 
condition may be more stringent that the positivity require- 
ment of (d 2 n/dA 2 )^ h . 

To present our results in a more concrete fashion, we study 
the phase diagram in different planes of the three-dimensional 
space spanned by the parameters (T, p, 1/ kpa). In particular, 
we have shown how the phase boundaries defined by the vari- 
ous stability conditions evolve with temperature. At relatively 
low polarization p < 0.2, there exists a stable superfluid phase 
at finite T, which may not be stable at T = 0. In a related 
fashion, we have shown where on the phase diagram pairing 
without condensation (that is, a pseudogap phase) appears. 

A general observation associated with these population im- 
balanced theories (which exclude the FFLO phase) is that in 
the important regime near unitarity, superfluidity only appears 
at intermediate temperatures 12411 . This was seen earlier in 
mean field approaches 13811 via studies of the quantity we de- 
fine as T^ IF , which was found to be double valued. In the 
present framework (which goes beyond by including pairing 
fluctuation effects), this intermediate temperature superfluid- 
ity appears via the presence of two T' c s at unitarity. Some- 
what above unitarity, and towards the BEC regime, we find 
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another mechanism for arriving at intermediate temperature 
superfluidity. This can occur due to a low temperature in- 
stability towards phase separation, associated with a negative 

signin(d 2 r!/dA 2 k h . 

While (d 2 il/dA 2 ) f2 ^h < is argued to lead to phase sepa- 
ration, one may imagine that a negative pair mass would result 
in a lower pair energy at a finite momentum than at q = 
and thus a possible LOFF-like condensate. On the other hand, 
a negative n s has also been argued |21] to be associated with 
a LOFF phase. Future studies are needed to unravel these var- 



ious possible phases within the unstable regimes in the phase 
diagrams. 
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